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1 Introduction 



Sum-frequency and difference-frequency generations are important processes 
in nonlinear optics (see, e.g., Ref. [TJ [2]). The aim of the present article is 
to develop a formalism describing such nonlinear optical processes in terms 
of Hertz potentials. The first part of the article is devoted to a derivation 
of the equations describing the generation of higher frequencies in a uniaxial 
nonlinear (up to second order) crystal in the paraxial approximation. Par- 
ticular attention is paid to the fact that, for extraordinary waves, the wave 
vectors do not coincide with the group velocity vector. This usually origi- 
nates some difficulties in optical experiments and must be taken into account 
carefully [3] . In the second part of the article, a possible scheme is proposed 
for generating higher frequencies using precisely this side-walk propagation. 
The idea is to combine phase matching, which involves wave vectors, with 
group velocity vectors. More specifically, the scheme consists in combining 
both ordinary and extraordinary beams in such a way that all three waves 
involved in the process are aligned in the same direction. The conditions 
to be fulfilled by the crystal parameters for this particular configuration are 
given explicitly. An analytic solution is also obtained for particular values of 
the amplitudes of the initial waves. 

The organization of this article is as follows. In section 2, a general for- 
mulation of the problem in terms of Hertz potentials is worked out following 
Nisbet's original treatment [4J. The results are applied in section 3 to the 
problem of sum-frequency generation. The possible alignment of the three 
group velocities involved in the scheme is studied in section 4. The evolution 
equations are presented in section 5, together with a particular analytic solu- 
tion. Finally, the particular case of second harmonic generation is considered 
in section 6. 



2 Propagation in a birefringent medium. 

Consider an anisotropic medium described by electric and magnetic field 
vectors, E and B, and electric displacement vector D. The Maxwell equations 
in the absence of free charges and currents (with magnetic permeability \i — 1 
and setting c = 1) are 

<9B 

V - B = , V x E + — = , (2.1) 
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<9D 

V-D = , VxB- — = 0. (2.2) 

dt y ' 

The effect of the material medium can be described by a polarization 
vector P such that 

D = E + 4vrP, 
the linear part being given in terms of the dyad 

?= ej_l + Ae s s, 

where s is the optical axis of symmetry of the medium, and e± and eii = 
e± + Ae are the permeability perpendicular and parallel to this symmetry 
axis respectively Then 

D(w, r) = e(co) E(w, r) + 4irP NL (u, r) 

where P^ is the non-linear contribution to the polarization vector. Here 
and in the following, Fourier transforms with respect to time of all quantities 
will be used. 

Following Nisbet[4j, the electromagnetic field can be described by two 
scalar Hertz potentials, to be called an d in the present paper, and 
two additional scalar potentials, to be called Uo and Ue- These potentials 
satisfy the equations 



V±U E + s x Wo = 47tP_ l jvl , (2.3) 

, V • efw) • V V E + e\\ (co)u 2 ^e - s ■ VU E = -terP\? L , (2.4) 

€±{uj) e±{u" 



V 2 ^o + e±(co)uj 2 ^o + iooUo = . (2.5) 

In these formulas, Vj_ is the gradient operator in the plane perpendicular to 
s. Eq. O implies 

U E = 4vr V ± • P ± NL , (2.6) 

\7\U = 4vr s • (V x P^) , (2.7) 

which permits to decouple the potentials ^o an d Uo from ^ E and Ue- 

As shown in a previous article (Hacyan and Jauregui|5j), the advantage 
of this formulation is that ^o an d ^e correspond to the potentials for the 
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ordinary and extraordinary waves respectively. The electromagnetic field is 
given by 

E= - ius x Wo + — V(s • V$ E ) + u 2 ^ E s - — VU E (2.8) 

and 

B = V x [V x (* oS )] - iuV x O b s) . (2.9) 

From these formulas, the two fundamental modes can be identified: the 
ordinary wave with s • Eq — and the extraordinary wave with s • = 0. 
The case P NL = with Uo,e — corresponds to the linear limit considered 
in 0. 

The non-linear polarization vector is usually defined as 

Pa L (u, r) = J dui J du 2 8{uj-u x - uj 2 )xabc(^i, u 2 )E h {u u r)E c (u 2 , r) 

(2.10) 

in the quadratic approximation, where Xabd^iT^i) is the (Fourier trans- 
formed) second-order susceptibility tensor (assumed to be homogeneous). 

A particularly important case is the one in which there is a discrete set 
of well defined frequencies u iy such that 

E a (u,r) = Y,^-uj l )E^(r) . 



Then the basic equations take the form 

V_L£/f (r) + s x VUg\r) = 47rP?(u; i , u k , r 

and 
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(2.11) 



e ll (u> i )u>i+— r ^V-e(GJ i yV 
for extraordinary waves and 



e ± (w;W 2 + V 2 



*«(r)+-^4 s-VC/f (r) = 4 7 rP i f ) (^,u; fc ,r) 



6±(W<) 



(2.12) 
(2.13) 



for ordinary waves, where (setting Ui = Uj + Uk) 

P®{u)i = uij + co k , r) = XabM = coj + u k )E^\r)E^{r) (2.14) 

and 

P^(u, =ui- u k , r) = XaUoOj =Ui- u k )E^\r)Ei k >(r) . (2.15) 
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3 Sum- frequency generation 



Consider a typical problem of sum-frequency generation. Suppose an or- 
dinary and an extraordinary waves, of frequencies oji and u 2 respectively, 
combine inside the crystal to generate an extraordinary wave of frequency 
cj 3 = u>i + ou 2 . Let ipi be the Hertz potential corresponding to frequencies u>i, 
and Ui the associated auxiliary potentials. Accordingly the basic equations 
take the form: 

^i 2 e±K) + V 2 1 i(r) = -tuJrU^v) , (3.1) 



ojle^(u 2 )e\\{uj 2 ) + V-?(u; 2 )-V 



(3.2) 



= e\\(u) 2 ) s • VC/ 2 (r) - 4ne ± (u 2 )Pu(u 2 = oj 3 - ou 1 ,r) , 



W3 2 e±(w 3 )e||(w 3 ) + V-eM-V # 3 (r) 



and 



where 



ei|(w 3 ) s • W 3 (r) - 47re_L(w 3 )P|f(w 3 =u 1 + u 2 ,r) 



V^C/i = 4vr s • [V x P\_(ui = uj 3 - oj 2 , r) 



ViC/ 2 = 4tt Vj 

h 



V 2 : C/ 3 = 47T 



= w 3 - w 2 ,r) 
P 2 (o; 2 = w 3 - wi,r) 



•P 2 (a; 2 
• P 3 (^ 3 

Xa6c(^l 
Xa6c(w 2 
Xa6c(w 3 



w 3 -a;i,r) , 
wi+w 2 ,r) , 

a; 3 -a; 2 )^(r)[P 2 (r)]* 
a; 1 +a; 2 )P fe 1 (r)P 2 (r). 



(3.3) 

(3.4) 
(3.5) 
(3.6) 

(3.7) 



As a next step, let us assume that the potentials have the form 

^ = ^(r)e lk - r (3.8) 

where 



k 2 = e±(ui) ujj 



(3.9) 
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and 

k i • • kj = e ± (u j )\s? j± + e\\(ujj)k\ = e ± (u j )e\\(u j ) uj) , (3.10) 
for j = 2, 3, and also 

\J X = Ul (r)e i(ks - k2) - p , (3.11) 

U 2 = u 2 (r)e l ^' k ^ r , (3.12) 

U 3 = M 3 (r)e i(kl+k2) - r . (3.13) 

In the above equations, A( r ) an d Wj(r) are slowly varying functions of r. 
Within this same approximation: 

Ei ~wi s x kii4i(r)e ikl " p (3.14) 

and 

W ?s - -4-t(s ■ k,)kJ^(r) e *V' - — ^rVC/, . (3.15) 



E ~ 



for j = 1,2. The last term in the above equation is quadratic in the elec- 
tromagnetic field; to be consistent, it must be neglected when evaluating the 
polarization vector up to second order. Accordingly: 

P 1 (u 1 =u 3 -u 2 ,r) = Pl A 3 (r)A5(r)e i < k3 - k2 >- r , (3.16) 
P 2 (w 2 = w 3 -wi,r) = p 2 A 3 (r)^(r)e^- kl )- r , 
P 3 (o; 3 = Wl +u; 2 ,r) = p 3 A^A^e^^ , 



where the vectors p^ are given in terms of Xa6c> s an d k, as 

Pl,a = Xabc(ui = ^3 - w 2 ) e 3jb e 2 , c , (3.17) 

£>2,a = Xabc{^2 = ^3 ~ ^i) e 3)b ei )C , 
£>3,a = Xabci^Z = ^1 + w 2 ) ei )6 e 2 , c , 



with 



for j = 2,3. 



ei = uis x ki , (3.18) 
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e j = t--t( s • kj)kj , (3.19) 

e ± {ujj) 
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The basic equations ( 13.1 p to ( 13. 6 p now take the form 



V 2 A l + 2ik x • VA X 



iAkr 



(3.20) 



V • e(w 2 ) • VA 2 + 2ik 2 ■ e(w 2 ) • VA 2 
«e||(a; a ) s • (k 3 - ki)w 2 - 47re ± (a; 2 )p 2 || A* X A 3 



-iAkr 



(3.21) 



V • e(w 3 ) • VA 3 + 2ik 3 • e(w 3 ) • VA 3 
ie||(w 3 ) s • (ki + k 2 )w 3 - 47re_ L (w 3 )p 3 ||A 1 A 2 



iAkr 



where Ak = k x + k 2 — k 3 , and 



and 



- (k 3± - k 2± )V = Am [s x (k 3 - k 2 )] • p 1± A* 2 A 3 

- (k 3± - k 1± ) 2 u 2 = Am (k 3± - k 1± ) ■ p 2± A\A 3 , 



- (ki ± + k 2± ) 2 M 3 = Am (ki ± + k 2 j_) • p 3± A X A 2 



(3.22) 



(3.23) 
(3.24) 

(3.25) 

within the same approximation (that is, keeping only terms of order k 2 ). 
3.1 Phase matching 

The phase matching condition is k 3 = ki + k 2 , in which case the above 
equations somewhat simplify: 



V 2 A 1 + 2zki • VA X 



-lUlUi, 



(3.26) 



V-e(u 2 )-VA 2 +2ik 2 -e(uj 2 )-VA 2 = ie\\(u 2 ) s-k 2 u 2 -47re ± (w 2 )p 2 ||^A 3 , (3.27) 
V-e(w 3 )-VA 3 +2ik 3 -e(w 3 )-VA 3 = ie\\(w 3 ) s-k 3 M 3 -4vre ± (w 3 )p 3 ||A 1 A 2 , (3.28) 
and 



- (ki ± ) V = Am (s x ki) ■ ViA* 2 A z 
-{k 2A _) 2 u 2 = Am k 2 ■ p 2 ±A{A 3 , 
-(k 3± ) 2 w 3 = Am k 3 • p 3 ±AiA 2 . 



(3.29) 
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Now, if no absorption is present, we have the following general relations: 

Xafec(^3 = Wi + U 2 ) = Xbca^x = U 3 - U 2 ) = Xcab(^2 = W 3 - Wi) 

(Kleinman [6]). Eqs. ( 13 . 1 7[) then imply 

ei ■ pi = e 2 • p 2 = e 3 • p* 3 = -C . (3.30) 

With this last condition, it follows after some lengthy but straightforward 
algebra [taking relations (I3.10p into account] that the basic equations (13.261 
EQgj) take the form 

An 

V 2 A l + 2*k x • VA 1 = - — - CA* 2 A 3 , (3.31) 

l k i-U 

An 

V • e(w 2 ) • VA 2 + 2zk 2 ■ e(w 2 ) ■ VA 2 = ^ e||(w 2 )e(w 2 ) CA*A 3 , (3.32) 

|k 2 ±| 

V ■ e(w 3 ) ■ VA 3 + 2zk 3 • e(u 3 ) ■ VA 3 = e,, (u 3 )e(u; 3 ) C*A 1 A 2 . (3.33) 

K3± 



4 Side-walk alignment 

Suppose we want to generate a frequency u 3 = wi + w 2 . Since rays 2 and 3 are 
extraordinary, they do not propagate along k 2 and k 3 , but rather along the 
directions e(2) • k 2 and e(3) • k 3 respectively, due to the side-walk effect (see 
e.g., Ref. [7]) . It it then possible to choose the directions of propagations 
in such a way that the three rays propagates along the same direction, in 
addition to the phase matching condition. This can be achieved setting 

ki + k 2 = k 3 , 

e(2) ■ k 2 = e ± (2)k 2± + e N (2)A; 2 ||S = ak x 

e(3) ■ k 3 = e ± (3)k 3± + e|,(2)fc 3 ||S = /3k x (4.1) 

(in this section, we set e(wj) — > e(i) in order to lighten the notation). This 
system of linear equations admits non-trivial solutions if the proportionality 
constants a and /3 take the values: 

en(2)e ± (2)Ae(3) 
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(4 . 2) 

where 

L> = e ± (3)Ae(2)-e ± (2)Ae(3). (4.3) 

Furthermore, it is evident that the three ray vectors kj and the optical axis 
s must be in the same plane. 

Dividing the wave vectors into components perpendicular and parallel to 
s, we have additionally the conditions 

fc?± + *i||=eL(l)wi 2 , (4-4) 



'' l -4^ 4- (4.6) 



ci|(3) e ± (3) 



and since 



h2L=-^k 1± , fell = -^y*i|| 



fcai=^* LL , fen = -^Jfe 1 „ 



it follows from ( 14. 5 p that 



Ae(2)fc? ± = (|^|y) 2 - e±(l)e±(2)«? , (4.7) 
Ae(2)^, = _(^.) 2 + e J .(i) C| |(2) W ? . (4.8) 



From (14.61) it also follows that 



/ tu 3 D \ ^ 



M3)^ = ( ^) j -e±(l)e±(3K, (4.9) 

Ae(3)^ = -(^|)) 2 + ^(1)^I|(3K • (4-10) 
Accordingly, the following relation is necessary for consistency: 



besides, of course, oj\ + u 2 = uj 3 . 

From the above formulas, it follows with some lengthy but straightfor- 
ward algebra that the angle 9\ between ki and s are given by the following 
equivalent formulas: 



sin 2 0i 



k 1± 



1 



/ uj 2 D y 
e ± (l)Ae(2) WiAe(3)/ 
1 / cu 3 D n2 



(4.12) 



e±(2) 

fcilP e ± (l)Ae(2) VwiAe(3)y Ae(2) 

^(3) 

e ± (l)Ae(3) \co 1 Ae(2)J Ae(3) 
[o; 2 Ae(2)] 2 e ± (3) - [a; 3 A6(3)] 2 6 ± (2) 
Ae(2)Ae(3) [wf Ae(3) - w 2 2 Ae(2)] ' 

The consistency conditions for these equations (since 1 > sin 2 6 > 0) are 



u; 3 Ae(3) en (2) > (<) u; 2 Ae(2) en (3) 



(4.13) 



if w 2 Ae(3) > (<)w 2 Ae(2). 



5 Evolution equations 



Eqs. (I3.311l3.33l) simplify considerably under the assumption that the only 
relevant spatial variations are along ki. Choosing the z axis along that 
direction, it follows that 

^-A 2 = C 2 A{A 3 , 
az 

^-A 3 = - CtA x A 2 , (5.1) 



dz 

where, using (14.1 j) and (14. 2 j) 



D 3 

° 2 = e 2 (2)Ae3(3) Cl (5 ' 3) 
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D 3 

° 3 = eJ(3)A e 3(2) Cl • (5 ' 4) 

Notice that all three coefficients C\, C 2 and C3 are complex, but have the 
same phase. 

From the above equations it follows that there is a conserved quantity: 
^(c 2 Ct\A x \ 2 + CzCl\A 2 \ 2 -C x C*\Az\ 2 ) = 0, (5.5) 

which is the Manley-Rowe relation [8]. 
A particular solution of Eqs. (15. ip is 

A\ = a\e lS sech(7^) 

A 2 = a 2 e lS sech(72) 

^3 = a 3e lS tanh(7^) , 

where 5 is the common phase of Cj and eij are real coefficients. Since the 
amplitudes of A\ and A 2 are given as initial conditions, the amplitude A3 of 
the generated wave follows from the relation 

a 3 = -\C 3 \a 2 a 2 ^ 1 , (5.6) 

with 

7 2 = |Ci||C r 3|ai=|C 2 ||C 3 |a?. (5.7) 

Thus the additional condition | C x I / 1 C2 1 = a i/ a 2 must be fulfilled for the 
above analytic solution to be valid. The ratio |Ci|/|C 2 | follows directly from 
Eqs. O and (Oil . 



6 Second harmonic generation 

Let us consider as a further example the generation of second harmonics by 
non-linear effects. Usually, under appropriate conditions, an ordinary wave of 
frequency u gives rise to an extraordinary wave of frequency 2u. Accordingly, 
the process is described by the equations given above, with the following 
identification: and k 2 correspond to the ordinary and extraordinary rays 
respectively, both with frequency u, and k 3 corresponds to the extraordinary 
wave of frequency 2u, that is: u\ = u 2 = u and 0J3 = 2oj. 
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In order to further lighten the notation, let us redefine e±(u>) = e and 
e_i_(2w) = e, and similarly for Ae and e\\. 

Then, according to the consistency condition (14. lip : 

from where it follows, using the definition of D, that 

e 



1 - 

e 



Ar "1 2 



2Ae 



(6.2) 



and therefore 



It also follows from ( 16 .2p that 



Ae 

< — < 4 . (6.4) 
Ae v ; 



This inequality must be satisfied in order to have triple alignment of the 
velocity vectors. 
Also 

. 2 Ae 2 e-4eAe 2 

sm &l = - — — —— t— r. (6.5) 

AeAe(4Ae - Ae) v ' 

Thus, if the optical axis s makes an angle <fi with the unit normal vector to 
the surface of the crystal, then according to Snel's law, 

sint = y / esin(6 l 1 — 0), (6.6) 

where i is the incidence angle to which the impinging ray must be directed 
in order to have a phase-matching assisted by side- walk alignment. Equation 
( 16 .4p must be satisfied. 

The evolution of the field is given by Eqs. (15. ip with its coefficient given 

by 

«- 

e 3 / Ae \ 6 

C ' = ij(2Ai) C " < 6 ' 7 > 
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7 Concluding remarks 



The formalism presented in this paper can be applied to other processes, 
such as difference-frequency generation and parametric down-conversion (to 
be considered in a forthcoming publication). As for the particular scheme of 
side-walk alignment herein proposed, it is left as a proposal to find crystals 
with the appropriate parameters, and to check its validity experimentally. 
Work supported by PAPIIT-UNAM project IN101511. 
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